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Abstract Wc discuss the possibility to obtain an elec- 
tromagnetic emission accompanying the gravitational 
waves emitted in the coalescence of a compact binary 
system. Motivated by the existence of black hole config- 
urations with open magnetic field lines along the rota- 
tion axis, we consider a magnetic dipole in the system, 
the evolution of which leads to (i) electromagnetic radi- 
ation, and (ii) a contribution to the gravitational radi- 
ation, the luminosity of both being evaluated. Starting 
from the observations on magnetars. we impose upper 
limits for both the electromagnetic emission and the 
contribution of the magnetic dipole to the gravitational 
wave emission. Adopting this model for the evolution 
of neutron star binaries leading to short gamma ray 
bursts, we compare the correction originated by the 
electromagnetic field to the gravitational waves emis- 
sion, finding that they are comparable for particular 
values of the magnetic field and of the orbital radius 
of the binary system. Finally we calculate the electro- 
magnetic and gravitational wave energy outputs which 
result comparable for some values of magnetic field and 
radius. 
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1 Introduction 

According to General Relativity, compact binary sys- 
tems (composed of neutron stars (NS), white dwarfs 
or black holes (BH)) are strong emitters of gravita- 
tional waves (GWs), the ripples of space-time due to 
the presence of accelerated masses (coalescence) . As a 
consequence of the gradual inspiralling, the system loses 
energy , linear a.nd an g ular momentum (IHulse fc Tavlor 



(1975) 



Nice (200, 



Weisberg fc Tavlor 



Tavlor I 



Stairs 



(|2no4) 



(|l989t) : 

The coalescence of a 
compact binary system can be split into three dis- 
tinct, but not sharply delimited phases: the inspi- 
ral, the plunge/merger and the ring-down. For the 
inspiral phase an accurate analytical description via 
the so-called Post-Newtonian (FN) expansion has been 
developed (for a review see ( Blanchet ( 20061 ) V Spin 
and mass quadrupolar contributions to this dynam - 
ics are also well-known ( Barker fc O'Connelll (1975I)') . 
with inclusion o f leadi n g order spin-orbit |Rvan (|l996 ) ; 



Kidder et al.l 



1997); 



Willi 




Maiar fc VasthI (l200a):ICornish fc Shapiro Key! (120101): 



19951): lO'Connelll (I2004D : 



Wang & WilJI(l2UUy|l:IArun et al. (l2UUy| ) :IKiem Je 
20101 ): lApostolatos et al.l |l995l IiQQgI ) : lMaia I^ (1200' 



Rieth 



Zeng fc Willi (120071) 1. spin-spin 
( Wang fc Willi (l2007^:lA^n et al.l (l2009 l) :lKlein fc Jetzeij 



and mass qua drupole mass monopole (jPoisson (^1998{); 
Racinj |2008l )) couplings. 

Similarly to the inspiral, the ring-down phase admits 
a perturbative analytical model which describes the 
damped oscillations of the single c ompact obje c t resu lt- 



(200i)). 



ing from the binary coalescence (jBerti et al. 
The plunge/merger phase is however non-perturbative, 
having no analytical description for generic systems, 



Templates on binarv inspiral (Buonanno fc Damoui 


(1999 


2000): Buonanno et al. (20031): Buonanno fc Damourl 


(2006 


):ICutler et al.l(ll993l): Damour et al.l( 


1998,l2000l): 


Damour (2001): Damour et al. (2002, 2003!) 


Finn fc Chernof* 
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( 1993f ): iFlanagan fc Hughed |l993l ): |Pan et all |2004l) ) 

and robust search algorithms were de veloped for GWs 



(20091): SenEUDta ( 


2009|); 


di Serafino & Riccio 




ometers (Arun et al. 


(I2OO9II 



di Serafino et al.i ( 2QfO.) ; 



and space-based interf er- 
Lang fc Hughed (|2009t) V 



In addition to GWs e mission, these syste m s emi t 



electromagnetic r adiati on (jCapozziello et al. I (|201l[ ): 



Capozziello et alj ([20101)). This emission is due to dy- 



the individual magnetic momenta of a binary neutron 
star system, we will introduce a global measure of the 
T^cctromagnctic properties of the system, given by a 
single magnetic dipole. 

Gravitational radiation transports away energy, lin- 
ear and orbital angular momentum from the system. 
While the total angular momentum (defined as the sum 
of the orbital angular momentum and possible individ- 
ual spins of the NS) exhibits instantaneous changes 
during orbital evolution, its secular evolution over 
timescales larger than the orbital period is a pure 
decrease in magnitude, w ith the direction kept fixed 
( Apostolatos et al.l |l994l) ). 

Having in mind (i) the configuration of open field 
lines aligned with the spin of the final compact object 
resulting from the coalescence, and (ii) the property of 
the gravitational radiation to conserve the direction of 
1992 |) : ISinha fc Mukhopadhvavl (|2010[) : lKouveliotou et ahh ^^he total angular momentum on timescales longer than 

the orbital period, we chose the magnetic field aligned 
with the total angular momentum of the system. 

In Sec. 2, we present an overview on GWs emission 
from coalescing binary systems. In Sec. 3, wc introduce 
and discuss the magnetic field and we also calculate the 
correction induced in the GW luminosity and ampli- 
tude. For comparison, in Sec. 4, we present the lumi- 
nosity of the accompanying electromagnetic radiation. 
These are compared and the conclusions presented in 
Sec. 5. 



namics of system and to behavior of plasma around 
the system. There are many different systems that 
emit gravitational and electromagnetic waves, and in 
many case the electromagnetic emission is very ener- 
getic (AGNs, QSOs, GRBs). 

Pulsars with surface magnetic fields in the range 
Q magnetars with surface ma gnetic fields up 
to 10^^^^^ G have been reported (.Duncan fc Thompson 



19981): [Thompson fc DuncanI (1X9951): lUsovl (19921) : 



Vasisht fc Gottheld (Il997l): 



Woods et al 




summarized in (jSinha fc Mukhopadhvavi (|201( 
upper limit of 10^^ G was also advanced from the re- 
quirement to avoid exotic s tructu re in the neutron stars 
( Sinha fc Mukhopadhvav ( 2010l )). An interaction of 
the magnetic dipoles modifies the gravitational radia- 
tion leaving the system, and hence the back-reaction 
to the orbit is slightly different. Limits on magnetar 
magnetic fields render these corrections to the second 
or even higher PN order. Neutron star binary coa- 
lescence, leading to a final black hole with a torus is 
a configuration believed to be the progenitor of short 
GRBs. In the black hole-torus central engine part of 
the torus matter can be conv erted into radiati o n, lea d- 



ing to the observed bursts ([Belczynski et al.l (|200 



Fully general-relativistic simulations of binary neutron 
star coalescen ce leading to th is confi guration have been 
presented in ( Rezzolla et al. ( 2010l )). with a detailed 
and accurate study of the torus. It was found that un- 
equal mass neutron star binaries led to a larger torus 
and an empirical formula in terms of the mass ratio 
and total mass was advanced for the torus mass. Nev- 
ertheless, the effects of the mag netic fields were not i n- 
cluded in the treatment of Ref. ( Rezzolla et al. ( 2010f )). 
Many efforts arc done in order to study the electro- 
magnetic emission from these systems, in particular 
from a computational point of view. It was showed 
how the dynamics of the binary systems induce the 
electr omagnetic emission and i ts lin k with GWs emis- 
sion JPalenzuela et al. I |2009l I2OI0I ): [Anderson et al. 

imm). 

In this paper we investigate a simple model for the in- 
spiral, which also takes into account the magnetic field. 
Rather then monitoring the complicated evolution of 



2 Gravitational v^^ave emission from coalescing 
binaries 

Let us summarize some basic concepts s of the gravi- 
tational radiation generated by the evolution of a com- 
pact binary system. A wide discussio n about the GW s 
emission can be found in literature (Maggior i (l2007l) 



([M , 

Misner et all (|l973l) : IWeinbergI ([197^). The starting 



point for any discussion of GWs is the approximation 
to General Relativity (GR), when the space-time met- 
ric can be treated as dev i ating only slightly froni a fla t 
metric iMaggiorel t00l\ ): IShapiro fc Teukolskvl |l983l ): 



Misner et al 



|l973 )). We start from linearized Einstein 



equation 



Dh„ 



167rG 



(1) 



where □ = rj^yd^^d'^ is the d'Alembertian, expressed in 
terms of the trace reversed metric perturbation 



1 



(2) 



where h^^ — g^^ ~ rj^^ and h = rj^^h^'^ (such that h ~ 
—h). Due to the imposed Lorentz (harmonic) gauge 



3 



dyh =0, only 6 of the 10 components of are 
independent. By setting T^^ = we obtain the vacuum 
wave equation: 



0. 



(3) 



GWs thus propagate at the speed of light. 

We evaluate the leading order contribution to the 
metric perturbations under the assumption that the in- 
ternal motions of the source are slow compared to the 
speed of light. We also assume that the source's self- 
gravity is negligible. Let us first introduce the momenta 
of the mass density 



M 



1 J d3xr0«(i,x), 



(4) 
(5) 
(6) 



and impose the conservation law d^jT^'^ = 0, valid in 
linearized gravity. For sources having strong gravita- 
tional field, as NSs or BHs, the mass density depends 
also on the binding energy, however for weak fields and 
small velocities T°°/c^ reduces to the rest-mass density 
p. Denoting the field h^^, which satisfies th e trans verse 



p. i^CllUUlllg UllV. /(-^i/ VV11H.11 OdUlJllCO Ull O UIOIIO VCIOC 

and traceless gauge conditions (jWeinbergl (jl972[ )). as 



transvcrse-traceless tensor hj^ , it is convenient to com- 
pute from it explic itly the two indipe ndent polarization 
states /i+ and /ix (jBuonannd (|2007l )). 

In terms of the traceless quadrupole tensor 



(7) 



{M. denoting the trace of Mij), the leading order ex- 
pression of the total radiated power is 



P 



-(Q 



(8) 



Henceforth, we will apply this formalism for a compact 
binary system with masses mi and m2, total mass AI = 
mi -f- m2 and reduced mass /i = mim2/(mi -|- m2). 
With the plane of motion chosen as the x-y plane and 
in the circular motion approximation the reduced mass 
particle has the coordinates r — (x, y, z) given by 



x{t) = r cos Lo t , y{t) ^ r sincot , z{t) ~ , 



(9) 



r being the relative distance between the two bodies. 
With p = p6{r) we find AP^ = px"^ x^ with the nontriv- 



ial components 

Mil = i ^r^ (1 -f cos2a;i) 

1 2 / 

M22 ~ (1 — cos2a;t) 

1 2 

From which we obtain 



4G 22(1 + cos^ ' 



cos{2ujt) 



AG_ 



p up' cos Q sin(2a; t) 



(10) 



(11) 



(12) 



Here the time i s shifted in order to eliminate the de- 
pendence on (/) ( Buonann"ol ( 200?! )'). From Eq. ([8]) the 
total radiated power emerges as, 



P = 



32 G p^ 



(13) 



In order to compensate for the loss of energy by GW 
emission, the radial separation r between the two bodies 
must decrease. The orbital frequency and consequently 
the GW frequency also changes in time, and can be 
derived using from the balance equation 



■bit 



dt 



-P. 



(14) 



3 Gravitomagnetic corrections on gravitational 
wave emission 

In this section we introduce the magnetic field with the 
desired structure, namely aligned with the total angular 
momentum J of the system. 

3.1 Gravitomagnetic-electromagnetic analogy 

The gravitomagnetic potential characterizing the com- 
pact binary is — (0, A), with 



J X r 



a- 



(15) 



where a = Gl c iTartaglia k. Ruggierol (2002), and r is 
the separation between the bodies of the binary system. 
Due to the fact that the motion of the binary system 
is entirely confined into the orbital plane, we evaluate 
all the physical quantities of interest in the plane it- 
self. This procedure is analogue to that used in GWs 
emission calculations, where one reduces the problem 



4 



only to the planar relative motion. The gravitomag- 
netic field is found as 



B = V X A = 



(16) 



We have used the assumption, that neutron stars ro- 
tate slowly, therefore the proper spin contributions to 
J are negligible, also that the relativistic (PN and 2PN 
corrections) are aligned to the Newtonian orbital angu- 
lar momentum, thus J • r = (for a motion confined to 
the x-y plane the only nonvanishing component of J is 
= J). Similarly, a gravitoelectric field emerges as 



E = —V X J 



(17) 



and they combine to a gravito-electromagnetie tensor 



F„ 



da An — dp An 



The remarkable property of this gravitomagnetic 
field B is its alignement to J, a property we seek 
for the magnetic field in our model. Therefore we 
will introduce an electromagnetic field analogous to the 
gravitoelectric and gravitomagnetic fields, by changing 
a = G/c into a = /io/47r, where fiQ is the vacuum mag- 
netic permeability. In what follows, by a we mean this 
expression and we refer to A, B, E, i^^^ as true electro- 
magnetic quantities. 

3.2 Gravitational wave energy loss and polarizations 

We characterize the source by the total stress-energy 
tensor 



^ ~ ^PF ^ ^EM 



composed by a perfect fluid part 



and an electromagnetic part 



(18) 



(19) 



(20) 



In order to calculate the momenta ([I])-® we need the 
00-components of these tensors. The 00-component of 
the perfect fluid contribution in a comoving frame de- 
flncd as 



u" = (u°, 0,0,0) , uqu° = 
simply becomes 

^00 _ rpPF _ 2 



(22) 



while for the electromagnetic contribution in the system 
in which the magnetic field is along the z-axis [such that 
J = (0, 0, J)] we find 



rpOO 

''■EM 



f-EM — 



J 



2c 



2 «J 2 2 



(23) 
(24) 



Despite apparences, T^j^j does not depend on masses, 
but it was written this way for an easy PN order eval- 
uation: 0{J/^ircf = 0{v/cf = 0{GM/c'^r). We 
also remark, that Oiriojcf' = Oivjcf'^ therefore the 



respective term of the expression 
dropped. 
Inserting 



T 



PF 



T 



EM 



pc 



fi^a^esM 



could be safely 



(25) 



into ^ , the corrections to the gravitational waves emis- 
sion can be derived. In a binary system, the pres- 
ence of an additional dipole-electromagnetic field has 
to be taken into account. The correction can be de- 
rived by considering also the electromagnetic stress- 
energy tensor beside the perfect fluid one. stress energy 
tensor. Thanks to the gravito-magnetic ormalism dis- 
cussed above, we can consider the electromagnetic con- 
tribution under the same standard of the gravitational 
one: this means to take into account the 00-component 
of the total stress-energy tensor given by the gravita- 
tional part (the mass density) and the electromagnetic 
contribution. In order to get the momenta (obtained 
by integrating the total stress-energy tensor in the vol- 
ume where the binary system is contained) we have to 
compute an integral consisting of a gravitational and 
an electromagnetic part: 



Mij = / (Pxp (x) XiXj + 



Gravitational part 



"^^W"^ - 2c ) x,xj 



Electromagnetic part 



(26) 



where p is the reduced mass of the binary system. We 
obtain 



(21) M, 



M 



1 + cos {2ujt) sin {2ujt) 
sin {2ut) 1 - cos (2a;i) 




(27) 



where is its trace, given by 



2c' 



'12 

-r ix) 



(28) 
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The nontrivial quadrupole tensor components are 
Qii = lM[l+3cos{2ujt)] , 

D 

Q22 = -M[l-3cos{2ojt)] , 



Q 



33 



n2 



-M sin (2wt) 



(29) 



Inserting them into Eq. (|14p we obtain the loss of en- 
ergy 



dE 
5c5 



dE 
H 



GWb 



(30) 



where 



stands to indicate the total energy 



GWt 



output of gravitational waves emission due to the pure 
gravitational term ) > ^nd the electromagnetic 



dE 



GW 



correction ( — — 1 . The polarizations result to be 

"-^ ' GWb 



1 



c^d \' 3c 
X [3 + cos(26')]cos(2a;i), 



(31) 



3c" 



^EM ) cos6'sin(2a;i). 



(32) 



Calculating h-^- and hx for a NS binary system using 

r = 10^ ^ 10^2^^, d-lO^^m, mi=m2 = l,4M0, 

after a few calculations we obtain that the corrective 
electromagnetic term is initially negligible with respect 
to the gravitational one, but it becomes not negligible 
when we increase the orbital radius, as well as the mag- 
netic field, as shown in Fig [T] Similar results can be 
obtained also for an equal-mass BH or White Dwarfs 
binary system for which the corrective terms are neg- 
ligible. However, it is possible to obtain for a binary 
system with a very strong dipole magnetic momentum, 
an electromagnetic contribution to GWs emission which 
is comparable to the standard gravitational one. 



4 Electromagnetic emission by the magnetized 
binary 

Now wc calculate the electromagnetic energy emission 
by binary systems in order to compare it to GWs energy 
emission. We consider the electromagnetic tensor com- 
ponents in order to calculate the time energy emission: 



-2a±{t) 



aJ 

cr^ 



Tf = -2ay (i) , 3 



aJ 
cr 



and finally 



Tf = 0. 



(33) 



(34) 



(35) 



Calculating the Poynting vector P, whose components 
are ~ cT^, it is possible to obtain the electromag- 
netic energy emitted in the unit-time as 



dE 



EM 



dt 



P-dS. 



where is the surface of the reduced mass object, 
which gives 



dE 



EM 



dt 



= -2acruj( ] S,, 



(36) 



From eq. ([T5)) . it is possible to get the dipole magnetic 
momentum J: 



J 



Br^ 
~2^' 



(37) 



Now we are able to compare the different contributions 
to energy emission: the standard GWs term, the elec- 
tromagnetic correction to it and then the pure electro- 
magnetic one. The results are shown in Tab. [TJ 



5 Conclusions 

GWs science has entered a new era and the recent years 
have been characterized by several major advances. 
For what concerns the most promising GW sources for 
ground-based and space-based detectors, notably, bi- 
nary systems composed of NS and BHs, our under- 
standing of the two-body problem and the GW gen- 
eration problem has improved significantly. The best- 
developed analytic approximation in General Relativ- 
ity is undoubtable the post-Newtonian method. In this 
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Table 1 The columns of the table, from left to right, indicate the radial separation between the two bodies (r), the orbital 
frequency {uj), the electromagnetic energy emitted in the unit-time (^~^~~^ ' ^'-'^^ '^^ energy in gravitational emission 

and its electromagnetic correction • "^^^ ranges are for magnetic fields B = (10^^ -^ 10^*) G. 



r (m) 


UJ (Hz) 




V dt J EM 


(crg/s) 






_ (dE) 

\ dt ) GWb 


(crg/s) 


10« 


3.04 X 10 


-4 


1.14 X 10^1 ^ 


1.14 X 10^3 


6.72 X 10^9 


6.8 


8 X 10^2 ^6. 


88 X 10^4 




9.61 X 10 


-6 


3.60 X 10^" ^ 


3.60 X 10^2 


6.72 X 10^4 


6. 


B8 X 10^ ^6. 


B8 X 10^1 


10" 


3.04 X 10 


-7 


1.14 X 10^0 H- 


1.14 X 10^2 


6.72 X 10^^ 


6. 


B8 X 10^ ^ 6. 


B8 X 10^*^ 


10^2 


9.61 X 10 


-9 


3.60 X 10^^ ^ 


3.60 X 10^1 


6.72 X 10^^ 


6. 


B8 X 10^ ^ 1.64 X 10^5 



paper, we have calculated the electromagnetic correc- 
tions to the GWs emitted by a coalescing binary system. 
We have considered an electromagnetic dipole-type field 
and calculated the electromagnetic contribution to the 
stress-energy tensor. We have obtained a correction to 
the standard gravitational energy-loss which becomes 
null setting the magnetic field to zero. In general, the 
electromagnetic correction term is negligible with re- 
spect to the gravitational one. This result holds also for 
black hole and white dwarfs binary systems. However, 
there could be coalescing binary systems with a large 
dipole magnetic momentum and a large orbital radius 
where the electromagnetic correction to GWs emission 
is relevant. In particular, the electromagnetic energy 
emitted by these binary systems is comparable to the 
gravitational wave output for given values of the mag- 
netic field and of the orbital radius. Revealing such a 
phenomenon by observations could be an electromag- 
netic signature for the gravitational waves emission. 




(c) 
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